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Abstract 

We study the phase transition from two different topological phases to the ferromagnetic phase 
by focusing on points of the phase transition. To this end, we present a detailed mapping from such 
models to the Ising model in a transverse field. Such a mapping is derived by re-writing the initial 
Hamiltonian in a new basis so that the final model in such a basis has a well-known approximated 
phase transition point. Specifically, we consider the toric codes and the color codes on some various 
lattices with Ising perturbation. Our results provide a useful table to compare the robustness of the 
topological codes and to explicitly show that the robustness of the topological codes depends on 
triangulation of their underlying lattices. 


1 Introduction 

One of the important problems in the control of quantum systems is the decoherence which emerges 
from undesirable interactions m. Specially in quantum computation processes, it is very important 
to find a quantum memory which is robust against emergent errors. In recent years, finding a solution 
for the decoherence problem has been one of the biggest challenges mia. One of the options to 
circumvent this problem is topological quantum computation where topological nature of a physical 
system is used to perform a fault tolerant quantum computation ala. 

The key property in the topological quantum computational model is topological order. Topological 
order is a new phase of matter which for the first time was presented by Wen in 1980 0. It is com¬ 
pletely different from ordinary order such as ferromagnetic order so that the topological orders cannot 
be described by Landau paradigm and symmetry breaking theory 0. Also there is no local order 
parameter which recognizes topological phase. In fact topological order is described by a non-local 
order parameter and it is independent from local details of physical systems 10. There are many 
physical systems in condensed matter physics with topological order 0 , for example in quantum Hall 
effects mni , high-temperature superconductors HU [El El and frustrated magnetism |E1[I5][I61. 
Among physical systems with topological order, lattice models have been taken into more consider¬ 
ation, because of their simplicity. In 1997, Kitaev presented a model Hamiltonian on a torus which 
could encode the subspace of two qubits in it’s ground states 0. Thereafter, topological color codes 
were introduced by Bombin et al. ifTTl [TSl [191 l20l where the unitary Clifford gates could be topologi¬ 
cally applied while this was not possible in the toric code. 
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Degeneracy of the ground state in the topological codes depend on topological nature of the model 
and local small perturbations cannot lift such a degeneracy. The robustness of the degeneracy of the 
ground state is only guaranteed in small perturbations so that when perturbation increases, at a critical 
value, a phase transition occurs. In fact, there is a topological phase transition (TPT) where the point of 
phase transition is a measure of the robustness of the code. The robustness of the topological codes are 
studied against thermal perturbations and local perturbations 112111^ 1^ l24l 1^ |2^ . As an example, 
robustness of the toric codes against parallel and transverse magnetic helds have been studied where 
magnetic held can be treated as an extrinsic factor which perturbs topological order ll27l l28l l2^ [30l . 
If the perturbation is as Ising interaction, there is a phase transition between topological phase and 
the ferromagnetic phase. Such a problem has been studied for the toric code in an Ising perturbation 
(TCI) on a square lattice ED and for the color code in Ising perturbation (CCI) on the honeycomb 
lattice in EH . 

In this paper, we focus on the effect of the Ising perturbation on topological codes on various lat¬ 
tices. We show that such models are mapped to Ising model in transverse held (IT) by re-writing the 
Hamiltonian in a new basis. Since the approximated phase transition points for the IT on different 
lattices are well-known EHEHESl . we provide a useful table to compare the phase transition points 
of topological codes on various lattices. 

We consider dependence of the phase transition point and also the robustness of topological codes on 
triangulation of their underlying lattices. In fact, for a pure topological model on a specihc lattice, 
there are some topological properties which are independent of triangulation of the lattice a. When 
the perturbation is added, the system is not topological. However, before the phase transition occurs, 
topological properties are still independent of triangulation. Since the point of phase transition is a 
property of the perturbed topological code, it is expected that the phase transition point depends on 
triangulation. We derive points of phase transition for perturbed topological codes on various triangu¬ 
lations and explicitly show dependence of phase transition points on triangulation of the lattice. 

In the first step, we compare the robustness of CCTs on two specific lattices. We show that they 
are mapped to Ising model in transverse field on two different lattices with different phase transition 
points. The same result is derived for the TCTs on three different lattices. We emphasize that although 
it is well known that the phase transition point in an ordinary quantum phase transition depends on the 
underlying lattice, such a problem has not been explicitly studied for a TPT. 

In the second step, we compare the robustness of the TCI and CCI. Such a work has been performed 
for the topological codes in a parallel field. In ll24ll . the authors have compared the robustness of the 
color code and the toric code in parallel field. They have concluded that the color code is more robust 
than the toric code against parallel magnetic field. We show that underlying lattices of these codes are 
important in such a problem so that the robustness of TCI on square lattice is less than the robustness 
of the CCI on honeycomb lattice while the TCI on honeycomb lattice is as robust as the CCI on hon¬ 
eycomb lattice. So the robustness of topological codes is very sensitive to triangulation of the lattices 
and we cannot compare different topological codes without considering the underlying lattices. The 
structure of this paper is as follows; 

In Section (|2]) we review topological properties of two topological codes, the color codes and the toric 
codes. In Section © we apply Ising perturbation to the topological color codes on various lattices. 
Then we apply the same perturbation to the toric code. In this section we present a detailed mapping 
to re-write Hamiltonian of the models in a new basis so that the initial models are converted to the 
Ising model in the transverse field. Finally, in Section (IHi, we provide a table of the phase transition 
points and compare the robustness of the topological codes. 

2 Topological code states 

In this section, we review some important properties of the topological quantum codes. Specifically, 
we consider two well-known topological codes, the color codes and the toric codes. Topological order 
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of these models is a non-local property which does not depend on the underlying lattices. For example 
the color codes can be defined on the honeycomb and square-octagonal lattices so that, in both cases, 
topological order is the same. 

2.1 Topological color codes 

In this subsection, we review the topological color codes (TCC). It has been shown that they can be 
used for topological implementation of Clifford unitary gates and it is an advantage of these codes 
compared with the toric codes US). 

The TCC models are defined on three-colorable lattices where one can color all faces (plaquettes) 
of the lattice by three different colors so that any two neighbor plaquettes are not in the same color 
m. Here we consider the TCC model on a honeycomb lattice, but we emphasize that topological 
properties are the same for all TCC models. To this end, consider a honeycomb lattice on a torus where 
spins live on vertices of the lattice, see Figure O- There are two kinds of operators corresponding to 
each plaquette of the lattice as follows: 

P. = X{X, , = ( 1 ) 

ieP iGP 

where X and Z denote the Pauli operators and az, respectively. If the number of the plaquettes is 
N, the number of the plaquette operators will be equal to 2N. The most important problem is to find 
the stabilized subspace of the code which is defined as follows: 

c = m \ PM = \i,),Pz\i^) = m. ( 2 ) 

This subspace can also be considered as the subspace of ground states of a Hamiltonian as: 

H = (3) 

p p 

By the fact that Px{i + Px) = (1 + Tk), it is simple to check that the following state is a ground state 
of the Hamiltonian ([2i: 

i^)=n(i+^-)ioo-o)- (4) 

p 

where |0) is eigenstate of the Pauli operator Z. Also there are many constraints between the plaquette 
operators. As it has been shown in figure ([T]i, the honeycomb lattice is colored by three colors such 
as blue, green, red. According to the periodic boundry conditions on the torus, it is clear that the 
following constraints will be hold on the plaquette operators: 

n =n =n n =n =n 

PGP PGG PGP PGP PGG PGP 

where P G R refers to the plaquttes which are colored by red and so on. Hence, from 2N plaquette 
operators, only 2N — 4 number of them are independent. Because the number of qubits is also equal 
to 2N, degeneracy of the ground state is sixteen-fold. 

In the following, we show that all ground states are constructed by applying some non-local operators 
to state (0. For defining these non-local operators, we consider that all edges of the honeycomb lattice 
can also be colored by three colors so that each edge with a specific color connects two plaquettes with 
the same color, see figure([T]i. Then we define three kinds of color strings on the lattice where each 
string moves on the edges and plaquettes with the same color and two endpoints of the strings live in 
the plaquettes, see figure O. We denote these strings by S^, . We define two string operators 

corresponding to each string in the form of: 

= n = n (6) 
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Figure 1: (Color online) In the color code on the honeycomb lattice, spins live on vertices of the 
lattice. The lattice is colored by three different colors, red, blue and green. Each link connects two 
plaquettes with the same color and it is colored by the same color, so all links are colored by three 
colors, red, blue and green. A green string move on green links and it has two end points in two green 
plaquettes. 


where the Pauli operators are applied on all qubits which live on the string S'^ and C can be each 
one of three colors, red, green or blue. If two end points of a string are connected to each other, it 
generates a closed loop and we will have a corresponding loop operator. Since each plaquette shares in 
two spins with the loop operator, it is simple to show each loop operator commutes with all plaquette 
operators in the Hamiltonian Q. By attention to topology of the torus, the loop operators are divided 
to two general classes. One class contains trivial loops which are contractible and another contains 
non-trivial loops which are non-contractible. 

Important point is about non-trivial operators. There are six fundamental non-trivial loops which are 
denoted by three colors and two directions on the torus, see figure (|2|i. Related to each kind of the 
non-trivial loops, there are two operators which are constructed by Z or AT operators. We denote these 
operators as: 

= n = n (7) 

,tT ^^^C,CT 

where Ic^a refers to a non-trivial loop and C denotes color of the loop and a = {0,1} denotes 
direction of the loop around the torus. An important property of the non-trivial loop operators is that 
they cannot be written as product of the plaquette operators in the Hamiltonian Q. But it is simple to 
show from three loop operators and only two number of them are independent and 

one can be written as product of two other ones and some plaquette operators. So we select only two 
colors of these loops, for example red and green. 

Because of the topology of the torus, two non-trivial loops where one is in direction 0 and another is 
in direction 1 cross each other, see Figure (|2l), so there is an anti-commutation relation between the 
non-trivial loop operators as follows: 




( 8 ) 


Other relations between these operators are commutation. According to the above relations the non¬ 
trivial loop operators generate a 2'^ dimensional space which can encode the space of four qubits. Also 
we can explicitly construct all sixteen bases of this space by the non-local operators in relation In 
fact, it is clear that the state |'(/’) as © is stabilized by Lfbecause they commute 
with operators then they are applied to state |00...0) which is eigenstate of the Pauli operator Z. 
So all bases of the degenerate subspace are constructed as follows: 

(9) 

where i,j, k and I are 0 or 1. We should emphasize that degeneracy of the ground state in this model 
is related to the topology of torus and the colorable structure of the lattice. 
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Figure 2; (Color online) Schematically we show three colored loops which wind around torus in two 
different directions. 



Figure 3: (Color online) In the toric code on the honeycomb lattice, spins live on edge of the lattice, 
vertex operators are applied on three spins which are neighbors of each vertex and plaquette operators 
are applied on six spins which are on the edges of each plaquette. Two non-trivial loops on lattice has 
been shown, one moves on edges of lattice which has been denoted by red string and another moves 
on edges of dual lattice which has been denoted by blue string. 


2.2 The toric codes 

The most simple of the topological codes are the toric codes which can be defined on each oriented 
graph Here we consider a honeycomb lattice on a torus with N spins which live on the edges of 
the lattice. The toric code on such a lattice is defined by the following plaquette and vertex operators: 

Bp = Y[Z, , A,=Y[x,, ( 10 ) 

zGp iGv 

where i € p refers to spins which belong to plaquette p and i € v refers to spins which are neighbor 
of vertex v, see Figure The stabilized subspace of this code can treated as the ground state of a 
Hamiltonian as: 

Hk = -J2Bp-J2Av (11) 

p V 

Since Ay{l + Ay) = 1 + Ay, it is simple to show that the following un-normalized state is a ground 
state of this Hamiltonian: 

!</.) = n(l + 2l„)|00...0). (12) 

V 

Also there are two constraints between the plaquette and vertex operators in the Hamiltonian as 
Ay = I and Bp = I. So the number of independent operators is equal to N — 2 and the 
ground state has four-fold degeneracy. To find other ground states we should consider two non-trivial 
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loops Z*’, around two directions of the torus where loops can be defined on the edges of lattice as 
well as the edges of dual lattice where the vertices, edges and faces of dual lattice are in one to one 
correspondence with the faces, edges and vertices of the lattice, respectively. There are two operators 
corresponding to a non-trivial loop l'^ on the lattice and a non-trivial loop on the dual lattice, see 
figure (O, in the form of: 

L-'" =l[X, , L-'" = Y[ Z,. (13) 

On one hand these operators commute with the Hamiltonian so they have joint eigenstates, and fur¬ 
thermore there exist the following anti-commutation relations between these operators: 

= 0 = 0. (14) 

Other relations between these operators are commutation. According to the above relations all four 
ground states of the model are derived as: 

( 15 ) 

where i,j are 0 or 1, and \(j>) is a state as in Ea. (fT2l i. Thus in the toric codes, topological order of 
the ground states is associated to the periodic conditions of the torus and it is independent from the 
underlying lattice. 


2.3 The robustness of the topological codes 

The degenerate ground state of the topological codes can be used as a robust memory for quantum 
computation goals. For example, consider a perturbation as magnetic field on spins in the toric code. 
Such a perturbation leads to splitting between the degenerate states. The important point in this case is 
that two ground states of the model are only mapped to each other by a non-local operator as product 
of Pauli operators on a non-trivial loop which is related to the topology of the model. Therefore, if the 
perturbation is small enough, such a splitting cannot happen. 

It is clear that the robustness of the degeneracy is lifted when perturbation increases. In fact a pertur¬ 
bation even in the first order leads to a transition term between the ground state and the exited states 
in the perturbed energy levels. This problem leads to a phase transition from topological order to an 
ordinary order. Finally, the topological order is lifted at a phase transition point. Therefore, the point 
of phase transition is a measure of the robustness of the topological codes. 


3 Ising perturbation on the topological code states 

In this section, we consider the topological codes on different two-dimensional lattices with Ising per¬ 
turbation. We show that these models are mapped to the Ising model in transverse field by re-writing 
the initial Hamiltonian in a non-local basis. In a general form, the Hamiltonian of the topological code 
with Ising perturbation is in the form of: 

H = —Jhxop — khj (16) 

where hrop is the Hamiltonian of the topological code and hi is the Hamiltonian of the Ising perturba¬ 
tion. At the first, let us consider two special regimes of the couplings. In case that J = 0 or K = oo, 
we have the Ising model which have a two-fold degenerate ground state where Z components of all 
spins are up or down. In this case system is in the ferromagnetic phase and a local order parameter 
can describe phase of the model E]. Now we consider a case that J = cxd or iT = 0, so we have the 
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topological model where the system is in a topological phase and any local order parameter cannot 
describe such a phase 0. Since in these two different regimes we have two different phases, it is a 
reasonable expectation that at a critical ratio of k/ J a phase transition occurs. 

We re-write the Hamiltonian (fTST l in a new basis and we show that in this new basis such a Hamiltonian 
is mapped to the Ising model in the transverse field on a 2D lattice. 


3.1 The topological color code on the honeycomb lattice 

Our method in studying the topological codes in presence of the Ising perturbation is to re-write 
Hamiltonian of the model in a new basis. It leads to convert the initial model to an Ising model in the 
transverse field. Similar method has been already used in other papers about the toric code on square 
lattice El and the color code on the honeycomb lattice 1^ . In this subsection, we explain the main 
idea for the TCC model on the honeycomb lattice then we apply it to other models on various lattices. 
To this end, we consider a Hamiltonian as: 


= (17) 

P P {i,3) 

where (i, j) denotes neighbor spins on the honeycomb lattice. In this model, we have considered a 
positiv coupling constant J for the TCC Hamiltonian and a positiv coupling constant K for the Ising 
perturbation. 

At the first, by the fact that the Ising perturbation commutes with the operators in the TCC model, 
we find that ground state of the new model is in a subspace which is stabilized by all operators P^. 
We re-write the Hamiltonian ( fTTl i in such a subspace by choosing an un-normalized basis as follows: 

= 11(1 + (-!)’■-P,)|00...0). (18) 

P 

where N is the number of plaquettes of the lattice and Pp’s are binary numbers which are attributed to 
each plaquette. In fact we can attribute a virtual spin to each plaquette where Vp denotes the value of 
the virtual spin in a plaquette p. The states in EadTSll are non-local so that they correspond to anyonic 
excitations of the TCC model. 

The number of states as ( fTSl l is 2^ and all of them are stabilized by the operators P^. About the 
above basis, we should emphasize that all these states are not independent, in fact there are con¬ 
straints between the plaquette operators in different colors according to relation Q. For example 
product of all operators P^ and P^ is equal to the Identity and we find that P^-Px\4>ri,r2,r3,...,rN) = 
\4>ri,r2,r3,...,rN)- This basis would be complete if we add some new states which are constructed by 
applying the non-trivial loop operators on the state \4>). 

In the following we should re-write the Hamiltonian (fTTl i in the basis (fTSl l. As we explained, all states 
(fTSl l are stabilized by P^’s, so Pz\(l>ri,r 2 ,r 3 ,...,rN) = l'/’r-i.r 2 .r- 3 ,...,rjv) ^nd since the number of P^’s is 
equal to N, we find that I.rs. 7-3 ,...) = ^I</'ri ,r 2 . 7-3 ....)■ 

Then we apply the operator P^ on the states (fTST l. By the fact that Pa;(l + {—ly^Px) = (—l)’'p(l + 
{—ly^Px), ws Px\4'ri,r2,r3,...,rN) = (~ I'/’ri,r 2 .r- 3 .....rw)■ We Can describe this final re¬ 

sult in an interesting view so that effect of Px in the basis (fTSl l is equivalent to effect of the Pauli 
operator Zp in a basis as \ri,r 2 , ■■■,riy) on the virtual spins. 

Interesting situation occurs when we apply the Ising term on the basis (fTSl l. To this end, as we ex¬ 
plained in the previous section, we can color all links of the honeycomb lattice in three colors. Thus 
we divide the Ising interactions to three parts corresponding to the color of their links as: 

Pr= E + E E (19) 

{i,j)&G {i,j)eR 
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Figure 4: (Color online) Left; An Ising perturbation corresponding to a blue link does not commute 
with two plaqutte operators Px and P^ in two endpoints of that link. Right; The virtual spins which 
live in blue plaquettes are denoted by blue circles. An Ising interaction ZiZj between two real spins on 
blue links does not commute with two plaquette operators Px on two it’s end points and it is equivalent 
to an Ising interaction XpXp' between virtual spins which live on corresponding plaquettes. Finally 
we have an Ising model on a triangular lattice. 


where {i,j) £ B denotes sites i, j which are in two endpoints of a blue link and so on. For example, 
we apply the Ising interaction on the blue links to the state \4>ri,r2,ri,,...,rM)■ Each link (i, j) has four 
neighbor plaquettes and it is simple to show that corresponding operator to (z, j) as ZiZj does not 
commute with two plaquette operators Px corresponding to two plaquettes in two endpoints of that 
link, see FigurellH left). Therefore, we find that ZiZj{l + {—lY^Px) = (1 + {—lYp~^^Px)ZiZj, 
where Px is corresponding operator to one of the end points of the edge {i,j) ■ Finally, we conclude 
when the Ising interaction corresponding to a blue link of the lattice is applied to basis | (pn , 7-2 .rs,... .r-w) > 
it rises up the virtual spins rp, rpi corresponding to two blue plaquettes in two end points of that link. 
Therefore, we can describe the Ising interaction ZiZj as an Ising interaction XpXpi on the basis 
|ri, r 2 ,..., tn), where p,p' are two blue plaquettes which are connected by a blue link. As you can 
see in Figure (0] right), if we denote this Ising interaction between the neighbor virtual spins by a link, 
we will have a triangular lattice where all virtual spins in blue plaquettes live on the vertices of that 
lattice. 

The above argument can also be repeated for the Ising interactions on the green and red links where, 
corresponding to each color, we will have an Ising model on a triangular lattice. 

By considering all terms in the Hamiltonian (fTTl) in the new basis, we will have a new Hamiltonian in 
the form of; 

H = -JN- {JZp + KXpXp,)- Y [JZp + KXpXp,)- Y [JZp + KXpXp,), (20) 

{p}eB {p}eG 

where {P} £ B refers to all virtual spins on the blue plaquettes and P, P' refer to the neighbor 
plaquettes in the same color on a triangular lattice. For finding the ground state of the Hamiltonian, 
we use this fact that the final Hamiltonian in (l20l) on the virtual spins has been written as summa¬ 
tion of three independent Hamiltonian on three triangular lattices with different colors. Therefore, it 
is enough to find the ground state of each part, independently. Interesting result is that each part of 
the Hamiltonian is an Ising model in a transverse field on a triangular lattice which a quantum phase 
transition has been known for it at J/iT Ri 4.77 1331. At J K limit all virtual spins are in positive 
direction of Z and correspondingly the real spins are in ground state of the color code. At K ^ J 
limit All virtual spins are in the ground state of the Ising model on the triangular lattice and corre¬ 
spondingly the real spins are in the ground state of the Ising model on the honeycomb lattice. 
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Figure 5: (Color online) In the color code on Square-octagonal lattice, spins live on vertices of the 
lattice. Similar to the honeycomb lattice, we color all plaquettes and links by three different colors, 
red, blue and green. 


3.2 the color code on the square-octagonal lattice 

The color codes can be dehned on any three-colorable lattices. Dehnition of the stabilizers are the 
same for all of them so that there are two kinds of stabilizers corresponding to each plaquette , Pz 
which are as product of X or Z operators on the spins of each plaquette. 

In this subsection, we consider the color code on a square-octagonal lattice in a Ising perturbation, see 
Figure (|5]l. In order to hnd the phase transition point, similar to the previous section, we re-write the 
Hamiltonian of the model in the basis (flSl l. All of arguments are similar to the previous section so 
that here we divide the Hamiltonian to three parts corresponding to three colors of the plaquettes of 
the lattice. Each part of the Hamiltonian in the basis (fTsT l is as an Ising model in the transverse held 
on the virtual spins in plaquettes with the same color. Only distinction is kind of the Ising interaction 
pattern on the virtual spins in each color. 

As it has been shown in Figure (|6]l, links of the lattice with green color generate a square lattice where 
the virtual spins live on it’s vertices. Moreover, blue links generate a square lattice where the virtual 
spins live on it’s vertices, see Figure©. Also we should be care that there are two blue links between 
each two blue neighbor plaquettes. This hnal point leads to an important distinction so that two Ising 
interactions in original model on the real spins correspond to an Ising interaction in the square lattice 
on the virtual spins. By considering a similar argument for the red plaquettes, the hnal Hamiltonian 
on the virtual spins will be as follows; 

H = -KN- {JZp+KXpXp,)- Y {JZp + 2KXpXp,)- Y {JZp + 2KXpXp,), (21) 
{p}gG {p}&R {p}es 

where {P} € G refers to all virtual spins on the green plaquettes and P, P' are neighbor plaquettes in 
the same color which generate a square lattice. Each part of the above Hamiltonian is an Ising model 
in the transverse held on a square lattice. Eor the hrst part, phase transition occurs at J/K ss 3 041 
and for two other parts, phase transition occurs at JI2K « 3. 'When the Ising perturbation is increased 
by K, at K « J/6 for the red and blue sublattices a phase transition occurs from topological order 
to ferromagnetic order but the green sublattice stays in topological phase. At AT w J/3 topological 
order is completely removed and we have a ferromagnetic order. We can compare this result with that 
for the color code on the honeycomb lattice where phase transition occurs at K JjA.TT. 

3.3 The toric code on the honeycomb lattice 

In this subsection, we study the toric code on a honeycomb lattice with the Ising perturbation. To 
this end, we consider the nearest neighbors of each spin and we assume an Ising interaction between 
them. As it has been shown in Eigure ([8] left), if we illustrate each Ising interaction by a link, Ising 
interaction pattern is described by a triangle-hexagonal lattice which the spins live on it’s vertices. In 
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Figure 6: (Color online) The virtual spins which live in green plaquettes are denoted by green circles. 
An Ising interaction ZiZj between two real spins on green links does not commute with two plaquette 
operators Px on two it’s end points and it is equivalent to an Ising interaction XpXp> between virtual 
spins which live on corresponding plaquettes. Finally we have an Ising model on a green square 
lattice. 



Figure 7: (Color online) The virtual spins which live in blue plaquettes are denoted by blue circles. 
There are two Ising interactions ZiZj between real spins on two blue links which connect two blue 
plaquettes. By the fact that these two Ising interactions does not commute with two plaquette operators 
Px on two it’s end points, we find that they are equivalent to an Ising interaction XpXpi between 
virtual spins which live on corresponding plaquettes. Finally we have an Ising model on a square 
lattice. 


this new lattice the vertex operators of the toric code correspond to the triangular plaquettes and the 
plaquette operators of the toric code correspond to the hexagonal plaquettes. Finally the toric code 
with the Ising perturbation can be described by a Hamiltonian in the form of: 

77 =( 22 ) 

V p {i,j) 

By the fact that the operator ZiZj commute with the operators Bp and the non-trivial loop operators 
L<r,z ^ jjjg ground state of the above Hamiltonian is in a subspace which is stabilized by 

Bp’s and Therefore, we select the following basis to re-write the Hamiltonian (l22l i: 

+ (" 1)"” |00...0) . (23) 

V 

Then, we find new form of the Hamiltonian (l22l i in this basis. First, By the fact that Bp’s commute with 
all operators A^, we find that = ^\4‘ri,r2,...,rN)- Second, we apply the operators 

by the fact that A„(l + (-l)’’"A„) = (-l')’’”(l-|-(-l)’’”X), we find that = 

It is useful to describe the final relation in view of virtual spins. To this end, 
if we insert virtual spins on the vertices of the initial honeycomb lattice, the state \ 4 >ri,r 2 ,...,rN) 
treated as product state |ri, r 2 ,..., r^) on the virtual spins. Therefore, effect of the operator Ay on the 
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Figure 8: (Color online) Left; Ising interaction pattern has been shown by a triangle-hexagonal lattice. 
Right; We denote two vertex operators by two dark triangles. An Ising interaction ZiZj between two 
real spins does not commute with two vertex operators on two it’s endpoints and it is equivalent to an 
Ising interaction X^Xyi between corresponding virtual spins which are denoted by red circles. 


new basis is equivalent to effect of operator Zy on the virtual spin r„. 

Finally, we apply the Ising operator ZiZj to the basis (l23l l. To this end, we consider a special link on 
the triangle-hexagonal lattice. Since the corresponding Ising operator to such a link does not commute 
with two vertex operators in two it’s endpoints, see Figure dS] right), we find the following relation 
between each Ising operator and a vertex operator in one of the it’s endpoints; 

Z,Zj{l + i-iy^Ay) = (1 + {-lY^+^Ay)Z,Zj. (24) 

Therefore, effect of the operator ZiZj on the basis (|2^ leads to flip the virtual spins which live on 
two endpoints of the corresponding link. It means that the Ising interaction ZiZj on the real spins in 
the initial model is equivalent to an Ising interaction XyXyi on the virtual spins. 

In the following, we should And effect of all Ising interactions in the initial model. To this end, 
according to Figure (|9]l, we color vertices of the initial honeycomb lattice by two different colors, 
blue and green. By attention to the above argument, each Ising interaction between the real spins 
in the initial lattice, similar to what has been shown in Figure (0 right), is equivalent to an Ising 
interaction XyX^r between two virtual spins in the same color. Finally, the Ising interactions in the 
initial honeycomb lattice convert to the Ising interactions on two triangular lattices in two different 
colors and the Hamiltonian (l22l i in term of the virtual spins are wrote as follows; 

H = -JN-{J Zv + k XyXy,) - {J Zy + k (25) 

vGTb {v,v')€Tb {v,v')gTc 

where Tg (To) refers to the blue (green) triangular lattice. Therefore, the initial Hamiltonian is 
converted to two independent Ising models in the transverse held on triangular lattices, see Figure 
ill. As we pointed in the previous subsection, the phase transition point of this model occurs at 
ratio J/k ss 4.77. We can compare this result with the toric code on square lattice where the phase 
transition point is at J/fc Ri 6 liSTl . 


3.4 The toric code on the triangular lattice 

As a Anal case, we consider the toric code on the triangular lattice. In such a model, spins live on the 
edges of a triangular lattice and two kinds stabilizer operators are defined corresponding to plaquettes 
and vertices of the lattice, see Figure lfTol i. Ising perturbation is applied between neighboring spins 
where pattern of the Ising interactions has been shown in Figure lfTTl left). Similar to the previous 
subsection, we consider a non-local basis as; 

+ (-1)"” 7l„) |00.. .0) . (26) 
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Figure 9: (Color online) The virtual spins on vertices of lattice are denoted by two green and red 
colors. Ising pattern between real spins is converted to two independent Ising model on virtual spins. 



Figure 10; (Color online) Plaquette and vertex operators has been shown. 


In such a basis, the plaqutte operators Bp are equal to Identity and the vertex operators Ay are equal 
to Pauli operator Zy on virtual spins which live on vertices of the lattices. 

For Ising interactions in the initial Hamiltonian, as it has been shown in Figure lfTTl right), by the 
fact that each Ising operator ZiZj does not commute with two neighboring vertex operators Ay and 
Ay', we find that the Ising operator in the initial model is equal to an Ising operator XyXy' in the 
new basis. If we apply such a transformation for all Ising operators, see Figure lfTSl) . we will have a 
triangular lattice on virtual spins with Ising interaction XyXyi between neighbor spins. Also there 
is an important point that two Ising interactions in the initial model are mapped to an new Ising 
interaction on virtual spins. Finally, the initial model in the new basis would be an Ising model in a 
transverse field as; 

E = -kN - J^Zy -2k XyXy'. (27) 

{v,,v') 

The phase transition point for this model is at J/2k r; 4.77. 

4 Comments on results of the mapping 

A summary of results of the previous section has been shown in Table ([T]). According to our deriva¬ 
tions, two well known sets of topological codes on the most important lattices has been mapped to 
Ising model in transverse field on different lattices. 

First, we compare the robustness of the TCFs on different lattices. The TCI on square lattice is mapped 
to IT on square lattice while the TCI on the honeycomb and the triangular lattices are mapped to IT on 
triangular lattice with different phase transition points. Therefore, the robustness of the TCFs depends 
on triangulation of the underlying lattice. 

Second, we compare the CCTs on two different lattices. The CCI on honeycomb lattice is mapped to 
the IT on triangular lattice while the CCI on square-octagonal lattice is mapped to the IT on square 
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Figure 11; (Color online) Left: Ising interaction pattern has been shown by a triangle-hexagonal 
lattice. Right: Corresponding Ising perturbation to a special link does not commute with two vertex 
operators in two endpoints of that. 



Figure 12: (Color online) The toric code model in Ising perturbation is mapped to a Ising model in 
transverse field on virtual spins of a triangular lattice. 


lattice. Therefore, the robustness of the CCTs also depends on triangulation of the underlying lattice. 
Third, we compare the robustness of the TCTs with the CCTs. The robustness of the TCI on square 
lattice is less than the robustness of the CCI on honeycomb lattice while the robustness is the same for 
the TCI and CCI on honeycomb lattice. Therefore, triangulation of the underlying lattice is important 
in comparison of two different topological codes so that there is not any way to compare the robust¬ 
ness of two different topological codes independent of their underlying lattices. 


5 Discussion 

In this paper, we derived the points of phase transition for two important topological codes on various 
lattices with the Ising perturbation, the toric code and the color code. To this end, we presented a 
detailed mapping between such models and the Ising model in transverse field. Our results explicitly 
showed that points of the phase transition and also robustness of perturbed topological codes depend 
on triangulation of lattices. Also, by comparing the robustness of the toric code and the color code, we 
showed that there is no way to compare the robustness of two different topological codes independent 
of triangulation of the underlying lattices. 
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The topological code in Ising perturbation 

Ising model in transverse field 

K/J 

The color code 

on honeycomb lattice 
on square — octagonal lattice 

cm triangular lattice 
on square lattice 

« 0.209 

Ri 0.333 

The toriccode 

on square lattice 
on honeycomb lattice 
on triangular lattice 

on square lattice 
cm triangular lattice 
on triangular lattice 

Ri 0.166 

« 0.209 
a; 0.104 


Table 1: The corresponding Ising model to each Topological codes and the phase transition points 
have been shown. 
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